February 2, 2008 1:42 Proceedings Trim Size: 9in x 6in 



paper'archivc'01 



o 
o 

^ ; SUMMING LOGS OF THE VELOCITY IN NRQCD 

O ! AND TOP THRESHOLD PHYSICS * 

Q 

ANDRE H. HOANG 

Max-Planck-Institut fur Physik 
(Werner-Heisenberg-Institut) , 
Fdhringer Ring 6, 
80805 Miinchen, Germany 
' E-mail: ahoang@mppmu.mpg.de 

To achieve reliable predictions of the top-antitop threshold cross section at a future 
C 3 ' e + e — Linear Collider logarithms of the top velocity need to be resummcd. I review 

the issues that make this problem complicated and show how the task can be 
achieved by renormalization in an effective theory using the so called velocity 
renormalization group. The most recent NNLL order results are discussed. 
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1. Introduction 

The so-called "threshold scan" of the total cross section line-shape of top 
pair production constitutes a major part of the top quark physics program 
at a future e + e~ Linear Collider. 1 From the location of the rise of the cross 
section a precise measurement of the top quark mass with experimental 
uncertainties below 100 MeV will be gained, while from the shape and the 
normalization of the cross section one can extract the top quark Yukawa 
coupling y t (for a light Higgs), the top width or the strong coupling. 2 

In the threshold region, y/s ~ 2m t ± lOGeV, the top quarks move 
with nonrelativistic velocity in the cm. frame. Due to the large top width, 
r ( « 1.5 GeV, toponium resonances cannot form and the cross section is a 
smooth function of the cm. energy. This also means that non-perturbative 
effects can be neglected for predictions of the total cross section. It is 
convenient to define the top velocity by mtv 2 = y/s — 2m t . Because in 
the loop expansion one encounters terms proportional to [a s /v) n from the 
instantaneous exchange of n time-like gluons, one has to count v ~ a s 
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and to carry out an expansion in a s as well as in v. Schematically the 
perturbative expansion of the R-ratio for the total top-antitop threshold 
cross section, R — <Tti/°~n+n- > nas the form 

R = vJ2 (— ) x { 1 ( L °) ; «a, « (NLO) ; o?„ a s v, v 2 (NNLO) \ . (1) 

k=0 V ^ ' 

The expansion scheme in Eq. (1) is called fixed-order expansion, although 
it actually involves summations of the terms proportional to (a s /v) n to all 
orders. The scheme can be implemented systematically using the factoriza- 
tion properties of Non-relativistic QCD (NRQCD) established by Lepage 
etal. 3 The NNLO QCD corrections to the total cross section were calcu- 
lated already some time ago. 4 ' 5 ' 6 ' 7 ' 8 ' 9 ' 10 One of the important results was 
that in a threshold top quark mass scheme (the most common schemes 
in use are the kinetic, PS and IS mass schemes 11 ) the theoretical uncer- 
tainties for a top mass measurement are comparable to the experimental 
uncertainties. However, it was also found that the NNLO corrections to the 
normalization of the cross section were as sizeable as the NLO corrections 
leading to a theoretical normalization uncertainty of at least 20%. 13 It was 
shown 12 ' 13 that this normalization uncertainty does not affect significantly 
the top mass measurement since the latter is mainly sensitive to the energy 
were the cross section rises. However, this normalization uncertainty jeop- 
ardizes competitive measurements of the top width, strong top coupling, 
and the top Yukawa coupling. Figure la shows the vector-current-induced 
R-ratio cr(e+e" -> 7* -> tt)/a(e + e~ -> at LO, NLO and NNLO in 

the fixed-order expansion for typical choices of parameters and renormal- 
ization scales. a It is the discrepancy between the seemingly well behaved 
NLO curves and the NNLO results which is the worrying issue. 

One way to understand the large NNLO corrections is to recall that the 
ti pair at threshold is non-relativistic and that its dynamics is governed 
by vastly different energy scales, the top mass (m t ~ 175 GeV, "hard"), 
the top three-momentum (p ~ mv ~ 25 GeV, "soft") and the top kinetic 
or potential energy (E ~ mv 2 ~ 3 — 4 GeV, "ultrasoft" ) . The top width 
T t w 1.5 GeV, which protects against non-perturbative effects (and which 
is the main reason rendering the top threshold a precision observable) is of 
order E, but can be neglected for most of the following conceptual consider- 
ations. This hierarchy of scales is the basis of NRQCD factorization, 3 which 



a The contributions originating from Z-cxchangc are about an order of magnitude smaller 
and can be neglected as far as discussions on QCD uncertainties are concerned. 
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separates short-distance physics at the scale m from long-distance physics 
at the non-relativistic scales p and E. In its original formulation NRQCD 
has one quantum field for each of the non-relativistic quarks and antiquarks 
and one quantum field for the low-energy gluons. NRQCD matrix elements 
and Wilson coefficients therefore involve three types of logarithmic terms, 




These logarithms cannot be rendered small (or summed up) for any single 
choice of the renormalization scale /x. However, this can be severe: for 
example, a s (m t ) \n(m^/E 2 ) ~ 0.8 for fi = m t in the case of top quark pair 
production close to threshold. Obviously, a more sophisticated effective 
field theory (EFT) framework has to be devised to resolve this issue and to 
allow to resum all logarithmic terms. 

In this talk I discuss some of the conceptual aspects of an effective the- 
ory that is capable of systematically summing all logarithms of the velocity 
that can appear in the description of the non-relativistic quark-antiquark 
dynamics. The tool needed is called the velocity renormalization group 
(VRG) 14 and it is the basis of a renormalization group improved perturba- 
tive expansion. To be specific, I will concentrate mostly on an EFT known 
as vNRQCD. 14,15,16,17 At thc cnd j wiU comc back to thc top threshold 

cross section in order to review the present status. 
2. The Proper Effective Theory 

To begin with, it is useful to summarize the properties one desires from an 
effective theory of non-relativistic quark-antiquark pairs: 

(1) IR (on-shell) fluctuations, including IR divergences, in the full the- 
ory are reproduced by EFT fields. This avoids large logarithms in 
matching conditions. 

(2) There is a well-defined and systematic power counting scheme (in 
v). It can be uniquely identified which operators to account for what 
to compute at a given order. Moreover, all EFT loop integrations 
are governed by only one single scale. 

(3) There is a consistent renormalization prescription to treat UV di- 
vergences and to formulate anomalous dimension which, eventually, 
will allow to sum all "large" logarithmic terms. 

(4) All symmetries (spin, gauge, etc.) are implemented. This gives the 
most predictive power. 
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(5) The Lagrangian is formulated in the regulator-independent way. 

The original NRQCD 3 does in fact not offer properties 2 and 3. This is 
because there is one single gluon field describing soft and ultrasoft gluon 
effects, while the power counting for soft and ultrasoft gluon effects differs. 
For summing logarithms of the type shown above this is a severe obstacle. 
For example, soft gluons are participating in the binding of the quark pair 
while ultrasoft gluons are responsible for retardation effects such as the 
Lamb shift, and their interactions have to be multipole expanded. The 
situation is, however, even more complicated because, although v <C 1 
ensures that there is a strong hierarchy m 3> mv 3> mv 2 , the soft and 
ultrasoft scales turn out to be correlated by the energy-momentum relation 
of the massive top quarks, E ~ p 2 /m t . 

One approach in the literature to resolve these issues has been to ignore 
the correlation of soft and ultrasoft scales at the beginning and to account 
for the scale hierarchy m mv 3> mv 2 in the traditional Wilsonian (step- 
wise) way. One starts with the NRQCD by Lepage etal. having soft gluons 
and then "integrates out" soft effects at the scale mv. This results in non- 
local heavy quark interactions, the potentials, and interactions of the heavy 
quarks with the remaining ultrasoft gluons. The theory below mv is called 
pNRQCD 18 (p for "potential") and naturally avoids any double counting 
of gluon effects. In the NRQCD-pNRQCD approach there is no unique 
power counting. A v counting taking into account the energy-momentum 
relation of the massive quarks exists in pNRQCD, while in the intermedi- 
ate NRQCD theory a l/m counting is applied and the heavy quarks are 
static. A particular property of pNRQCD is that the correlation of the soft 
and ultrasoft scales has eventually to be implemented back into the the- 
ory when UV-divergences in massive quark loops occur. Renormalization 
group running caused by such divergences then also couples the ultrasoft 
renormalization scale together with the soft matching scale back to the hard 
scale m t . 19 

In an alternative approach one accounts for the correlation of the soft 
and the ultrasoft scales from the very beginning which then forbids using 
a step-wise approach for the hierarchy of the momentum and the energy 
scale. Rather, there is only one EFT below the scale m t which then con- 
tains simultaneously soft and ultrasoft gluons. This approach is called 
vNRQCD 14, 15,16,17 (v for "velocity") and allows for a consistent v power 
counting that obeys the correlation E ~ p 2 /m t at all scales below m t . 
There is an ultrasoft renormalization scale, nu, for loops dominated by the 
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ultrasoft energy scale, and a soft renormalization scale, jj,s, for loops dom- 
inated by the soft momentum scale which also includes the massive quark 
loops. Both renormalization scales are correlated, [ijj — fig/m = mv 2 . 
The running of coefficients and operators in the EFT is then expressed in 
terms of the dimension-zero scaling parameter v, where v = 1 corresponds 
to the hard matching scale and v = v is the low-energy scale where ma- 
trix elements are computed. The resulting renormalization group scaling, 
describing the correlated running coming from soft and ultrasoft effects is 
called the VRG 14 . As long as one does not have to account for UV diver- 
gences coming from massive quark loops there is a rather straightforward 
correspondence between the operator structure in vNRQCD and the one in 
the NRQCD-pNRQCD approach. 

The effective vNRQCD Lagrangian (defined in the cm. frame) is build 
from heavy potential quarks and antiquarks (ip p , x p ), soft gluons, ghosts, 
and massless quarks (A%, c q , ip q ) and ultrasoft gluons, ghosts, and mass- 
less quarks (A 11 , c, ip us ). Double counting for the gluon effects is avoided 
since ultrasoft gluons reproduce only the physical gluon poles where k° ~ 
k ~ mv 2 , while soft gluons only have poles with k° ~ k ~ mv. All soft 
loops are made infrared-finite and at the same time all ultrasoft divergences 
in ultrasoft loops are made correspond to the hard scale m by the pull-up 
mechanism. 20 It is essential that both soft and ultrasoft gluons are included 
at all scales below m because the heavy quark equation of motion corre- 
lates the soft and ultrasoft scales. The dependences on soft energies and 
momenta of the heavy quark and soft gluon fields appear as labels on the 
fields, while only the lowest-energy ultrasoft fluctuations are associated by 
an explicit coordinate dependence. Formally this is achieved by a phase re- 
definition for the potential and soft fields 14 , <f>{x) — ► J2k e ~ zk X( Pk{x), where 
k denotes momenta ~ mv and d^(j>k(x) ~ mv 2 (f)k(x). 



3. The Effective Theory Action 

The effective vNRQCD Lagrangian for a non-relativistic ti in its cm. frame 
in a angular momentum S-wave and color singlet state has terms 14 ' 15 ' 16 



iD° + i^-+dm t - 



^ P + O -> x) 



p k 



p,p',q,q,a 
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+ E \p" a p - p va v\ 2 + ■ ■ ■ E vfv(p, p') < VpxL^x-p 

p p,p' 

+ E 7^%/ ABC (P - p')-(^A c )[^ T ^ P X t _ p ^ B X-p] + • • •(?) 
p,p' lP PJ 

where color and spin indices have been suppressed and g s = g a (m t i l ), 
g u = g s (m t v 2 ), g s being the QCD coupling. All coefficients are functions of 
the renormalization parameter v, and all explicit soft momentum labels are 
summed over. Massless quarks and ghost terms are not displayed and there 
are other terms needed for renormalization purposes 21,22 not shown here. 
In d = 4 — 2e dimensions powers of and [i s are uniquely determined 
by the mass dimension and the v counting of each operator. The covariant 
derivative contains only the ultrasoft gluon field, = + in\ J g u A tl . The 
first line in Eq. (3) contains the heavy quark kinetic terms and their inter- 
action with ultrasoft gluons. There are 4-quark potential-like interactions 
of the form (k = (p — p')) 

vl i \ V c V k n 2 V r (p 2 + p' 2 ) v 2 , v s q2 

where S is the total ti spin operator. 

At NNLL order for the total cross section the coefficient V c of the 1/k 2 
potential has to be matched at two loops 20 because it contributes at the 
LL level, whereas the coefficients of the order 1/m 2 potentials only have to 
matched at the Born level 15 . The l/(m|k|)-type potentials are of order a 2 
and have to be matched at two loops. 21 There are also 4-quark interactions 
with the radiation of an ultrasoft gluon (last line) and interactions between 
quarks and soft gluons (second line). Here, due to momentum conservation 
at least two soft gluons are required. The sum of the potential terms shown 
in Eq. (4) and time-ordered products of soft interactions contribute to the 
instantaneous interactions between the top-antitop quark pair (that are 
traditionally called "potentials"). 

Particularly important for the phenomcnological application are the ef- 
fective operators involving the top quark width and the residual mass term 
Suit (first line). The width term arises by including the absorptive elec- 
troweak contribution from the W-b final state of the top 2-point function 
into the vNRQCD matching conditions at u = 1. The width is of order 
r t ~ E <~ D° <~ p 2 /m t ~ mi) 2 and has to be included in the heavy quark 
propagator. It is this width term which explicitly serves as an IR cutoff 
and which causes the smooth cross section shape visible in Figs. 1. Techni- 
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cally, it corresponds to a shift of the cm. energy into the positive complex 
plane by iT t ■ 23 Beyond LL order, however, the structure of electroweak 
corrections is more complicated. 9 ' 24 The residual mass term is non-zero in 
the threshold mass schemes, and cancels the order Aqcd ambiguity arising 
from the higher order QCD corrections to the Coulomb potential V c /k 2 
that would otherwise destabilize the perturbative behavior of the location 
where the threshold cross section rises. 13 (In the pole scheme Sm t — 0, and 
the destabilization due to an O(Aqcd) renormalon takes place. 13 ) 

Besides the interactions contained in the effective Lagrangian that de- 
scribe the dynamics of the tt pair we also need external currents that 
describe the production of the top quarks. For e+e~ annihilation these 
currents are induced by the exchange of a virtual photon or a Z boson. 
At NNLL order we need the vector S-wave currents Jp = c\{v)0 Pt \ + 

C2{v)O pa , where O pl = ^cr(ia 2 )X- p , O p2 = ^A^pM^)** p 
and the axial- vector P-wave current Jp = c 3 (v)O p .3, where O p ,3 = 
^ipp^[<T,cr ■ p](i<j2)x*- p - The currents O p2 and O p ^ lead to contribu- 
tions in the total cross section that are v 2 -suppressed with respect to those 
of the current O p \. Thus, at NNLL order, two-loop matching is needed 
for ci and Born level matching for c 2 and c 3 . 

4. The Computations 

The computations for the total cross section in the EFT are performed in 
the three steps, 

i) matching to QCD (including in principle also electroweak interac- 
tions) at \is = l^u = m (v = 1); 

ii) determination of the anomalous dimensions and running with the 
VRG from v = 1 to v = vo, were vo is the typical heavy quark 
velocity (v ~ 0.15 for tt threshold production); 

iii) computation of the EFT matrix elements at v = vq. 

Conceptually, the novel aspect of the computations is the summation of 
the logarithms of the correlated soft and ultrasoft scales using the VRG. 
Although these calculations are usual standard once renormalization is ex- 
pressed in terms of the renormalization parameter v, they are non-trivial 
conceptually. This is because the correlation of soft and ultrasoft renor- 
malization scales in the evolution is crucial to achieve the proper summa- 
tion of logarithms. This was demonstrated in QED, 25 where higher order 
logarithmic terms could only be reproduced correctly using the VRG. A 
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nice demonstration at the technical level on the difference between a cor- 
related and an uncorrelated running was given in 26 . Technically the ne- 
cessity of correlated running arises whenever divergences in massive quark 
loops and ultrasoft or soft loops contribute simultaneously to an anoma- 
lous dimension. It was also shown 17,22 that subdivergences in the NNLL 
(3-loop) renormalization of the current O p ,i, arising from diagrams where 
UV-divergences from massive quark and ultrasoft or soft loops arise simul- 
taneously in one single diagram, could be treated consistently in the VRG, 
while an uncorrelated treatment of soft and ultrasoft renormalization fails. 
Through steps i) and ii) all logarithmic terms involving v are summed into 
the operator coefficients and the matrix elements in iii) involve ln(v/v) 
terms that are not parametrically large. 

As far as QCD effects are concerned, for the renormalization group im- 
proved tt threshold cross section in vNRQCD all necessary components in 
the steps i), ii) and iii) are known at NNLL order with the exception of 
step ii) for the current coefficient c\. The running of the potential coeffi- 
cients Vi was determined in Refs. 15 > 20 > 27 > 21 and of the current coefficients 
d in Refs. 27 ' 21 . Up to notational differences there is agreement with the 
potential coefficients obtained in pNRQCD. 28 The running of c\ is fully 
known at NLL order, 27,21 while at NNLL order there is still some work to 
do. The full set of non-mixing contributions coming from UV-divergences 
in actual three-loop vertex diagrams have been determined in 22 . A non- 
trivial cross check for this computation, which checked the In v term in 
c\{v) was also carried out. 29 The full set of mixing contributions, which 
arise from the next-order running of the Wilson coefficients appearing in 
the NLL anomalous dimension of c\ is still unknown. Recently, the mixing 
contribution arising from the spin-dependent C(l/m 2 ) potential (last term 
in Eq. (4) was computed in NRQCD-pNRQCD. 30 



5. The Top Pair Total Cross Section at Threshold 

The total cross section for e+e~ — > 7*, Z* — > tt at threshold at NNLL order 
in renormalization group improved perturbation theory has the form 

( s ) — \_F V ( S ) K"( s ) + F a (s)R a (s)^ , (5) 

where F v,a are trivial functions of the electric charges and the isospin of 
the electron and the top quark and of the weak mixing angle. The vector 
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and axial- vector i?-ratios have the form 17 
4tt 



IV (a) 



Im 



s 

47T 



c 1 (i/) Ai(v, m, v) + 2 Ci(v) c 2 (y) Ai{v, m tl v) 



R a (s) = — Im c^(v) A 3 {v,m t ,iy) 



(6) 
(7) 



where the A^s are T-products of the effective theory currents described be- 
fore. As mentioned before, the running of c\ is only partly know at NNLL 
order. These T-products are related to the zero-distance Green function 




346 347 348 349 350 351 352 353 354 




346 347 348 349 350 351 352 353 354 
Vs(GeV) 

Figure 1. Panel a) shows the results for Qf K° with M 1S = 175 GeV and T t = 1.43 GeV 
in fixed-order perturbation theory at LO (dotted lines), NLO (dashed lines) and NNLO 
(solid lines). Panel b) shows the results for Q^R V with the same parameters in renor- 
malization group improved perturbation theory at LL (dotted lines), NLL (dashed lines) 
and NNLL (solid lines) order. For each order curves are plotted for v = 0.15, 0.20, and 
0.3. The effects of initial state radiation, beamstrahlung and the beam energy spread at 
a e+e~ collider are not included. The plots are from Ref. [32], 
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of the equation of motion of the 4-quark 4-point functions obtained in the 
EFT. (At NNLL order this is a common 2-body Schrodinger equation, since 
higher Fock quark-antiquark-gluon states only contribute to the renormal- 
ization of operators.) I use the IS mass 9,31 since it is the threshold mass 
that is most closely related to the peak position visible in the theoretical 
prediction. Other threshold masses are also viable 13 but not discussed here 
further. 

In Fig. lb I have displayed the photon induced cross section Q^R V at LL 
(dotted blue lines), NLL (dashed green lines) and NNLL (solid red lines) 
order for v = 0.15,0.2 and 0.3. Figure la shows the corresponding results 
in fixed-order perturbation theory already discussed earlier. Compared to 
the fixed-order results with the same scales the improvement is substan- 
tial, particularly around the peak position and for smaller energies, but 
the somewhat inconvenient behavior that the NLL and the NNLL order 
predictions do not overlap remains. At present the uncertainty of the nor- 
malization of the total cross section is 8a t t/a t i ~ ±6% due to this visibly 
large shift between the NLL and the NNLL order predictions. 32 In fact, an 
uncertainty not larger than 3% would be required to have the theoretical 
errors of measurements of the top Yukawa coupling, T t and a s match with 
the ones presently expected from experiments. 2 It is, however, premature 
to draw definite conclusions as long as the NNLL mixing contributions to 
ci have not yet been determined completely. 

6. Acknowledgement 

I would like to thank the organizers of the conference for the hospitality 
and the pleasant atmosphere at the Theoretical Physics Institute. 

References 

1. R. D. Heuer, D. J. Miller, F. Richard and P. M. Zerwas, hep-ph/0106315. 
T. Abe et al. [American Linear Collider Working Group Collaboration], in 
Proc. of the APS/DPF/DPB Summer Study on the Future of Particle Physics 
(Snowmass 2001) , ed. N. Graf, arXiv:hep-ex/0106057. 

2. M. Martinez and R. Miquel, Eur. Phys. J. C 27, 49 (2003) [arXiv:hep- 
ph/0207315]. 

3. W.E. Caswell and CP. Lepage, Phys. Lett. 167B, 437 (1986). G.T. Bodwin, 
E. Braaten and G.P Lepage, Phys. Rev. D51, 1125 (1995), ibid. D55, 5853 
(1997). 

4. A. H. Hoang and T. Teubner, Phys. Rev. D58, 114023 (1998) [hep- 
ph/9801397]. 



2, 2008 1:42 Proceedings Trim Size: 9in x 6in 



paper'archivc'01 



11 

5. K. Melnikov and A. Yclkhovsky, Nucl. Phys. B528, 59 (1998) [hep- 
ph/9802379]. 

6. O. Yakovlev, Phys. Lett. B457, 170 (1999) [hep-ph/9808463] . 

7. M. Beneke, A. Signer and V. A. Smirnov, Phys. Lett. B454, 137 (1999) 
[hep-ph/9903260]. 

8. T. Nagano, A. Ota and Y. Sumino, Phys. Rev. D60, 114014 (1999) [hep- 
ph/9903498]. 

9. A. H. Hoang and T. Teubner, Phys. Rev. D60, 114027 (1999) [hep- 
ph/9904468]. 

10. A. A. Penin and A. A. Pivovarov, Nucl. Phys. B550, 375 (1999) [hep- 
ph/9810496]; Phys. Atom. Nucl. 64, 275 (2001) [Yad. Fiz. 64, 323 (2001)] 
[hep-ph/9904278]. 

11. M. Battaglia et al, arXiv:hep-ph/0304132. 

12. D. Peralta, M. Martinez and R. Miquel, talk presented at the 4th Interna- 
tional Workshop on Linear Colliders, Sitges, Barcelona, Spain, April 28 - 
May 5 1999. 

13. A. H. Hoang et al, in Eur. Phys. J. direct C3, 1 (2000) [hep-ph/0001286]. 

14. M. Luke, A. Manohar and I. Rothstein, Phys. Rev. D61, 074025 (2000) 
[arXiv:hep-ph/9910209]. 

15. A.V. Manohar and I.W. Stewart, Phys. Rev. D 62, 014033 (2000) [arXiv:hep- 
ph/9912226]. 

16. A.V. Manohar and I.W. Stewart, Phys. Rev. D62, 074015 (2000) [arXiv:hep- 
ph/0003032]. 

17. A. H. Hoang, A. V. Manohar, I. W. Stewart and T. Teubner, Phys. Rev. D 
65, 014014 (2002) [arXiv:hep-ph/0107144]. 

18. N. Brambilla, A. Pineda, J. Soto and A. Vairo, Nucl. Phys. B 566, 275 (2000) 
[arXiv:hep-ph/9907240]. 

19. A. Pineda, Phys. Rev. D 66, 054022 (2002) [arXiv:hep-ph/0110216]. 

20. A.H. Hoang, A.V. Manohar and I.W. Stewart, Phys. Rev. D 64, 014033 
(2001) [arXiv:hep-ph/0102257]. 

21. A. H. Hoang and I. W. Stewart, Phys. Rev. D 67, 114020 (2003) [arXiv:hep- 
ph/0209340]. 

22. A. H. Hoang, Phys. Rev. D 69, 034009 (2004) [arXiv:hep-ph/0307376]. 

23. V. S. Fadin and V. A. Khoze, JETP Lett. 46, 525 (1987). 

24. A. H. Hoang and C Reifier, in preparation. 

25. A. V. Manohar and I. W. Stewart, Phys. Rev. Lett. 85, 2248 (2000) 
arXiv:hep-ph/00004018; 

26. A. V. Manohar, J. Soto and I. W. Stewart, Phys. Lett. B 486, 400 (2000) 
arXiv:hep-ph/0006096; 

27. A.V. Manohar and I.W. Stewart, Phys. Rev. D63, 054004 (2001). [arXiv:hep- 
ph/0003107]. 

28. A. Pineda, Phys. Rev. D 65, 074007 (2002) [arXiv:hep-ph/0109117]. 

29. B. A. Kniehl, A. A. Penin, M. Steinhauser and V. A. Smirnov, Phys. Rev. 
Lett. 90, 212001 (2003) [arXiv:hep-ph/0210161]. 

30. A. A. Penin, A. Pineda, V. A. Smirnov and M. Steinhauser, Nucl. Phys. B 
699, 183 (2004) [arXiv:hep-ph/0406175]. 



February 2, 2008 1:42 Proceedings Trim Size: 9in x 6in 



papcr'archivc'Ol 



12 



31. A. H. Hoang, Z. Ligeti and A. V. Manohar, Phys. Rev. Lett. 82, 277 (1999) 
[hep-ph/9809423]; Phys. Rev. D59, 074017 (1999) [hep-ph/9811239]. 

32. A. H. Hoang, Acta Phys. Polon. 34, 4491 (2003) [arXiv:hep-ph/0310301]. 



